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Abstract. The BBW method remains one of most demanded tool to derive full set of Cepheid 
astrophysical parameters. Surface brightness version of the BBW technique was preferentially 
used during last decades to calculate Cepheid radii and to improve PLC relations. Its implemen- 
tation requires a priory knowledge of Cepheid reddening value. We propose a new version of the 
Baade-Becker-Wesselink technique, which allows one to independently determine the colour ex- 
cess and the intrinsic colour of a radially pulsating star, in addition to its radius, luminosity, and 
distance. It is considered to be a generalization of the Balona light curve modelling approach. 
The method also allows calibration of the function F{CIo) = BC{CIo) + 10 x log{Tef f{CIo)) 
for the class of pulsating stars considered. We apply this technique to a number of classical 
Cepheids with very accurate light and radial- velocity curves. The new technique can also be 
applied to other pulsating variables, e.g. RR Lyraes. We also discuss the possible dependence of 
the projection factor on the pulsation phase. 

Keywords. Distance scale, Cepheid radii, reddenings and luminosities, period - luminosity 
relation, projection factor. 



1. Introduction 

Classical Cepheids are the key standard candles. They are used to set the zero point 
of the extragalactic distance scale (Freedman et al 2001) and also serve as important 
tracers of young populations (Binney and Merrifield 1998). They owe their popularity to 
their high luminosities and photometric variability (which make them easy to identify 
and observe even at large distances) and the fact that the luminosities, intrinsic colours, 
and ages of these stars are closely related to such an easy-to-determine quantity as the 
period of variability. 

It would be best to calibrate the Cepheid period-luminosity (PL), period-colour (PC), 
and period-luminosity-colour relations using distances based on trigonometric parallaxes. 
However, the most precisely measured parallaxes of even the nearest Cepeheids remain 
insufficiently accurate and, more importantly, they may be fraught with so far uncov- 
ered systematic errors. Here the Baade-Becker-Wesselink method (BBW; Baade 1926, 
Becker 1940, Wesselink 1946) comes in handy, because it allows the Cepheid distances 
(along with the physical parameters of these stars) to be inferred, thereby providing 
an independent check of results based on geometric methods (e.g., trigonometric and 
statistical parallaxes). The surface brightness technique was used most frequently and 
effectively during the past few decades. It is based on the relation between so-called 
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limb- darkened surface brightness parameter and the normal colours of Cepheids (Barnes 
and Evans 1976). Moreover, it critically depends on the adopted reddening value. 

Cepheid reddening values can be estimated from medium- and broad-band photo- 
metric observations (including multicolour PL relations) and from spectroscopic data 
(Dean, Warren and Cousins 1978, Fernie 1987, Fernie 1990, 1994, Fernie et al. 1995, 
Berdnikov, Vozyakova and Dambis 1996, Berdnikov, Dambis and Vozyakova 2000, An- 
drievsky et al. 2002a, b, Kovtyukh et al. 2008, Kim, Moon and Yushchenko 2011). AU 
these methods use proper calibrations and relationships between key stellar parameters. 
However, there exist large (up to 0.2mag) scatter in the estimates of the colour ex- 
cess for individual Cepheids, and the structure of the Cepheid instability strip is still 
vague. In the review on the Hubble constant and the Cepheid distance scale Madore 
and Freedman (1991) noted that "...any attempt to disentangle the effects of differen- 
tial reddening and true color deviations within the instability strip must rely first on 
a precise and thoroughly independent determination of the intrinsic structure of the 
period-luminosity-color relation." And next: "...independent reddenings and distances to 
individual calibrator Cepheids must be available." 

It should also be noted that reliable values of the colour excess and of the total- 
to-selective extinction ratio, say, Av/E(V — /), are extremely important when we use 
Wesenheit function to derive Cepheid luminosities from precise trigonometric parallaxes 
(Groenewegen and Oudmaijer 2000, Sandage et al. 2006, van Leewen et al. 2007). To 
convert Wesenheit index Wyi to the absolute magnitude. My, the intrinsic colour of the 
Cepheid, (V — /)o, and the proper value of Av/E(y — I) are needed. We suppose that 
self-consistent and independent reddening estimates for individual Cepheids can reduce 
underestimated systematical errors induced by large variations of the absorption law 
(Fitzpatrick and Massa 2007) and can even result in Ax law in optics and NIR. Therefore, 
the search for independent estimates of Cepheid reddening values is still actual, and this 
is our primary aim. 

Both BBW techniques - surface brightness (radius- variation modelling) and maximum 
likelihood (light-curve modelling) - are based on the same astrophysical background but 
make use somewhat different calibrations (limb-darkened surface brightness parameter, 
bolometric correction - effective temperature pair) on the normal colours. Here we pro- 
pose a generalization of the Balona (1977) light-curve modelling technique, which allows 
one to independently determine not only the star's distance and physical parameters, but 
also the amount of interstellar reddening, and even calibrate the dependence of a linear 
combination of the bolometric correction and effective temperature on intrinsic colour 
(Rastorguev and Dambis 2011). 



2. Theoretical background 

We now briefly outline the method. First, the bolometric luminosity of a star at any 
time instant is given by the following relation, which immediately follows from the Stefan- 
Boltzmann law: 

L/Lq = {R/Ref X {T/Tq)\ (2.1) 

Here L, and T are the star's current bolometric luminosity, radius, and effective 
temperature, respectively, and the subscript denotes the corresponding solar values. 
Given that the bolometric absolute magnitude, M^oZ, is related to bolometric luminosity 
as 



Mi,oi = MboiQ - 2.5 X log{L/LQ), 
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we can simply derive from Eq. (2.1): 

M^oi - Mi,oiQ = -bx logiR/Re) - 10 x log{T/TQ) (2.2) 

Now, Mi)oi can be written in terms of the absolute magnitude M in some photometric 
band and the corresponding bolometric correction, BC^ i.e. 

Mboi = M^BC, 

and the absolute magnitude M can be written as: 

M = m-A-bx log{d/10 pc). 

Here m, A, and d are the star's apparent magnitude and interstellar extinction in the 
corresponding photometric band, respectively, and d is the heliocentric distance of the 
star in pc. We can therefore rewrite Eq. (2.2) as follows: 

m = A + 5 X log{d/10 pc) + M^oiq + 10 x log{TQ) 

- 5 X log{R/RQ) - 5C - 10 X logiT). (2.3) 

Let us introduce the function F{CIq) = 5C + 10 x log(T), the apparent distance modulus 
(m — M)app = A-\-b X log{d/ 10 pc), and rewrite Eq. (2.3) as the light curve model: 

m = Y -5x log{R/RQ) - F. (2.4) 

where constant 

Y={m- M)app + MboiQ + 10 x log{TQ). 
We now recall that interstellar extinction, A, can be determined from the colour excess 
CEdisA = RxxCE, where Rx is the total-to-selective extinction ratio for the passband- 
colour pair considered, whereas Mi^oIq^ Rq, and Tq are rather precisely known quantities. 
The quantity F{CIq) — BC + 10 x log{T) is a function of intrinsic colour index CIq = 
CI — CE. Balona (1977) used a very crude approximation for the effective temperature 
and bolometric correction, reducing the right-hand of the light curve model (2.4) to the 
linear function of the observed colour, with the coefficients containing the colour excess 
in a latent form. 

The key point of our approach is that the values of function F are computed from 
the already available calibrations of the bolometric correction BC{CIo) and the effective 
temperature logT{CIo) (Flower 1996, Bessel, Castelli and Plez 1998, Alonso, Arribas and 
Martinez- Roger 1999, Sekiguchi and Fukugita 2000, Ramirez and Melendez 2005, Biazzo 
et al. 2007, Gonzalez Hernandez and Bonifacio 2009). These calibrations are expressed 
as high-order power series in the intrinsic colour: 

N 

F(C/o) = ao + ^afe•C/o^ (2.5) 

k=l 

with known {ak} and N ^ 7; in some cases the decomposition also includes the metal- 
licity, [Fe/H], and/or gravity {log g) terms. 

As for the stellar radius, R, its current value can be determined by integrating the 
star's radial- velocity curve over time, using dt = (P/27r) • dip: 

R{t) -Ro = -PF . r {Vr{t) - Vj) • (P/27r) • dip, 

J (po 

where Rq is the radius value at the phase ipo [we use mean radius, < R >= {Rmin + 
Rmax)/'^ ]; V-f the systemic radial velocity; ip the current phase of the radial velocity 



4 



A. S. Rastorguev et al. 



curve, P the star's pulsation period, and PF the projection factor that accounts for the 
difference between the pulsation and radial velocities. Given the observables (light curve 
and apparent magnitudes, m, colour curve and apparent colour indices, C/, and radial 
velocity curve and Vr) and known quantities for the Sun, we end up with the following 
unknowns: distance, mean radius, < R and colour excess, CE, which can be found 
simply using a maximum-likelihood technique (nonlinear optimization). 

For Cepheids with large amplitudes of light and colour curves {AC I ^ OAmag), it 
is also possible to apply a more general technique by setting the expansion coefficients 
{ak} in Eq. (2.5) free and treating them as unknowns. We expanded the function F = 
BC + 10 X log(T) in Eq. (2.4) into a power series of the intrinsic colour index CIq^ of 
a well-studied "standard" star (e.g., a Per or some other bright star) with accurately 
known T^^: 

N 

F = BC'^ + 10 X logiT'^) + ^ak-{CI-CE- CF^^f (2.6) 

k=i 

The best fit to the light curve is provided with the optimal expansion order N '^b — ^. 
We use this modification to calculate the physical parameters and reddening CE of the 
Cepheids, as well as the calibration F(C/o) = BC{CIo) + 10 x log{Teff{CIo)) for a star 
of given metallicity [Fe/H] and log g cite [(Rastorguev and Dambis 2011)]RD11. 

3. Observational data, constants, and best calibrations 

Cepheid photometric data were obtained by L. Berdnikov; his extensive multicolor 
photoelectric and CCD photometry of classical Cepheids is described in Berdnikov (1995). 
Very accurate radial- velocity measurements of 165 northern Cepheids were measured 
with Moscow CORAVEL spectrometer (Tokovinin 1987) during the period 1987-2011 
(Gorynya et al. 1992, 1996, 1998). The total number of individual measurements is about 
11,000 (the latest data are currently being prepared for publication). The photometric 
and spectral data sets are nearly synchronous, to prevent any systematic errors in the 
computed radii (up to 30%) and other parameters owing to evolutionary period changes 
resulting in phase shifts between light, colour and radial velocity variations. We adopt 
Tq = 5777 K, MboiQ = +4.76^ (Gray 2005). We proceeded from (V,B - V) data and 
found as the best solutions for the V-band light curve and {B — V) color curve those 
computed using the F[{B — V)o] function based on two calibrations (Flower 1996, Bessel, 
Castelli and Plez 1998) of similar slope; the poorer results obtained using the other cited 
calibrations can be explained by the fact that the latter involved an insufficient number 
of supergiant stars. 

4. The projection factor 

There is yet no consensus as to which projection factor (PF) should be used for Cepheid 
variables (Nardetto et al. 2004, Groenewegen 2007, Nardetto et al. 2007, Nardetto et 
al. 2009). Different approaches (constant or period-dependent PF values) lead to sys- 
tematic differences in the inferred Cepheid parameters, first and foremost in the radii. 

Bearing in mind the specific features of CORAVEL measurements, Rastorguev (2010) 
introduced phase-dependent PF. Its value is calculated from flux integration across the 
stellar limb and depends on limb darkening coefficient, e (linear darkening law D{(p) = 
1— e-he cosif, (f is the angle between line-of-sight direction and the surface element normal 
vector), and the photospheric velocity, dr/dt. The true line profile will be broadened by 
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Figure 1. Example: true line profile (e = 0.75, dr/dt = 40 km/s), observed line profile and its 
Gaussian approximation (with a spectrograph instrumental line width of 5*0 = 4 km/s). The 
measured Vr value is indicated by the arrow; the PF value is shown. 
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Figure 2. Projection factor, PF, as a function of \dr /dt\ for an instrumental width of 
aSo = 6 km/s and for different values of limb-darkening coefficient 0.5 < e < 1. 



any spectral instrument, and usually we approximate it by Gaussian curve to measure the 
radial velocity as the coordinate of the maximum (see Fig. 1). In this case, the measured 
radial velocity will additionally depend on the instrument's spectral linewidth, (so PF 
value should be "adjusted" to the spectrograph used for radial- velocity measurements), 
and on the photospheric velocity, \dr/dt\. 

PF values were estimated for 2.5 < \dr / dt\ < 60 km/s; 0<e<l; 4<6'o<8 km/s. 
We see, that the maximum of the normal approximation is shifted relative to the tip of 
the line profile. This shift depends on \dr/dt\ and So- Fig. 2 shows that the calculated PF 
values may differ considerably from the "standard" and widely used value PF = 1.31. 
Variations of PF have also been reported by, e.g. Nardetto et al. (2004). We provide useful 
analytic approximation for the projection factor, PF, as a three-parameter exponential 
expression, based on 10,000 numerical experiments (Rastorguev 2010): 

PF ^ ai • exp[—{dr/dt)'^/{2 ■ a2)] + as, where ai,a2,a3 are functions of e and Sq: 

ai ^ -0.068 • e - 0.0078 • 5'o + 0.217 

a2 ^ +1.69 • e + 2.477 • + 9.833 

as - -0.121 • e + 0.009 • So + 1.297 

Overall, the rms residual for this analytical expression is approximately 0.003. In prac- 
tice, the PF value for measured a given radial velocity should be determined by iteration 
for known values of So- To compare our results on Gepheids with other calculations, we 
finally adopted a moderate dependence of PF on the period advocated by Nardetto et 
al. (2007), although we repeated all calculations with other variants of PF dependence 
on the period and pulsation phase to assure the stability of the calculated colour excess. 
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Figure 3. Examples of the light-curve modelling. Crosses: observations; solid line: model. Fit 

quality: rms ^ 0.02 — 0.04ma^. 



5. Results and discussion 

To test the new method, we used the maximum-likelihood technique to solve Eq. (2.3) 
for the V-band light curve and B — V colour curve for several classical Cepheids residing 
in young open clusters: SZ Tau, CF Cas, U Sgr, DL Cas, and GY Sge, and found good 
agreement for the calculated reddening values with those determined for the host clus- 
ters (Rastorguev and Dambis 2011). A weak sensitivity of calculated reddening, Eb-Vi 
on the adopted PF value (constant or period/phase-dependent) is explained by the very 
strong dependence of the light curve's amplitude on effective temperature, ~ 10 x log(T)^ 
and - as a consequence - on the dereddened colour. Although the internal errors of the 
reddening Eb-v seem to be very small, the values determined using the two best calibra- 
tions (Flower 1996, Bessel, Castelh and Plez 1998), may differ by as much as 0.03-0.05^, 
due to the systematic shift between these two calibrations. Fig. 3 shows the final fit to 
the F-band light curves for some Cepheids of our sample. 

PL relation and the instability strip for our Cepheid sample are shown on Fig. 4. Note 
that the inferred radii and luminosities of large fraction of the Cepheids with P < 4^ are 
too large for fundamental tone pulsations; in most cases this may be indirectly evidenced 
by their low colour amplitudes. 

To refine the calibration F{CIq) (2.6), we tried to use a Per as the "standard" star, 
with T'^ ^ (6240±20) K, [Fe/H] ^ -0.28±0.06 (Lee et al. 2006), {B-Vy^ ^ 0.48^ and 
Eb-v ^ 0.09^ (WEBDA, for a Per cluster). To take into account the effect of metallicity 
on the zero-point F{CIoy\ we estimated the gradient dF{CIoy^ /d[Fe/ H] ^ +0.24 from 
the calibrations by Alonso et al. 1999, Sekiguchi & Fukugita 2000, Gonzalez Hernandez 
& Bonifacio 2009. In some cases (particularly for large- amplitude color variations) the 
"free" calibration, i.e. Eq. (2.6), can markedly improve the model fit to the observed 
light curve of the Cepheid variable. Fig. 5 shows the example of calibrations of the 
F functions derived from nine Cepheids with different [Fe/H] and log g values. The 
temperature scatter at Tg// ~ 6600 — 5100 K amounts to 3 — 5%. 

When applied to an extensive sample of Cepheid variables with homogeneous photo- 
metric data and detailed radial-velocity curves, the new method is expected to lead to a 
completely independent scale of Cepheid reddening values and to refine the PL relation. 
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Figure 4. PL relation (left) and the instability strip (right) for approximately 40 Milky Way 
Cepheids. Intensity mean absolute magnitudes and colours are shown. Most of short-period 
Cepheids seem to be overtone pulsators. Cepheids residing near the edges of the instability strip 
and labeled by their names all have very small amplitudes of the light and colour curves. 




Figure 5. Inferred calibrations of the F((B — V)o) function, Eq. (2.6), for nine Cepheids with 
large amplitudes of their colour variation. 
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